THE DOUBLE SHUFFLE RELATIONS FOR p-ADIC 
MULTIPLE ZETA VALUES 
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Abstract. We give a proof of double shuffle relations for p-adic 
multiple zeta values by developing higher dimensional version of 
tangential base points and discussing a relationship with two (and 
one) variable p-adic multiple poly logarithms. 
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0. Introduction 

In this paper we will prove a set of formulas, known as double shuffle 
relations, relating the p-adic multiple zeta values defined by the second 
named author in |Flj . These formulas are analogues of formulas for the 
usual (complex) multiple zeta values. These have a very simple proof 
which unfortunately does not translate directly to the p-adic world. 

Recall that the (complex) multiple zeta value C(k), where k stands 
for the multi- index k = (hi, ... , k m ), is defined by the formula 



(0.1) 



C(k) 



0<ni<-<n, 



n. 
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The series is easily seen to be convergent assuming that k m > 1. 

Multiple zeta values satisfy two types of so called shuffle product 
formulas, expressing a product of multiple zeta values as a linear com- 
bination of other such values. The first type of formulas are called 
series shuffle product formulas (sometimes called by harmonic product 
formulas). The simplest example is the relation 

(0.2) C(An) • C{k 2 ) = C(h, k 2 ) + C(h, h) + C(h + h) , 

which is easily obtained from the expression (j0.1|) by noting that the 
left hand side is a summation over an infinite square of pairs (721,712) 
of the summand in (jO.lJ) . and that summing over the lower triangle 
(respectively the upper triangle, respectively the diagonal) gives the 
three terms on the right hand side. Every series shuffle product formula 
has this type of proof. 

The second type of shuffle product formulas, known as iterated in- 
tegral shuffle product formulas, is somewhat harder to establish and 
follows from the description of multiple zeta values in terms of multiple 
polylogarithms. More precisely. The one variable multiple polyloga- 
rithm is defined by the formula 



(0-3) Lik(s)= E k~ ' 

0<ni<-<n m "-l lb m 

near z = 0. It can then be extended as a multi- valued function to 
P : (C) — {0, 1, oo}. We clearly have the relation lim^i Lik(z) = COO- 
Multiple polylogarithms can be written using the theory of iterated 
integrals due to Chen [CTj j. In other words, they satisfy a system 
of unipotent differential equations. This gives an integral expression 
for multiple polylogarithms. By substituting z = 1 and splitting the 
domain of integration in the right way we obtain the iterated integral 
shuffle product formulas, a simple example of which is the formula 
(0-4) 

ah)<(h) = £ ( h ah-i, k+o+x; ( kl ~ 1 + 3 )c(h-j, h 

i=0 ^ 1 ' 3=0 \ 3 / 

In |Flj the second named author defined the p-adic version of mul- 
tiple zeta values and studied some of their properties. The defining 
formula ()0.1|) can not be directly used p-adically because the defining 
series does not converge. Instead, one must use an indirect approach 
based on the theory of Coleman integration |Co| IBeslj . Coleman's 
theory defines p-adic analytic continuation for solutions of unipotent 
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differential equations "along Frobenius" . Coleman used his theory ini- 
tially to define p-adic polylogarithms. In |Flj Coleman integration was 
used to define one variable p-adic multiple polylogarithms. Taking the 
limit at 1 in the right way one obtains the definition of p-adic multiple 
zeta values. It is by no means trivial that the limit even exists or is 
independent of choices, and this is the main result of |Flj . 

Given their definition, it is not surprising that for p-adic multiple zeta 
values it is the iterated integral shuffle product formulas that are easier 
to obtain p-adically. In |Flj the series shuffle product formulas were not 
obtained. The purpose of this work is to prove fTheorem l7.1|) these for- 
mulas, and as a consequence the double shuffle relations (Corollary [73) 
for p-adic multiple zeta values. 

To prove the main theorem it is necessary to use the theory of Cole- 
man integration in several variables developed by the first named au- 
thor in F3esl a . The reason for this is quite simple - If one tries to 
replace multiple zeta values by multiple polylogarithms in the proof of 
(|U.2j) sketched above one easily establishes the formula 

(0.5) U kl (z) Li k2 (w) = hi klM (z,w) + Li k2>kl (w,z) + U kl+k . 2 (zw) , 

which is a two variable formula. It seems impossible to obtain a one 
variable version of the same formula. The proof of the main theorem 
thus consists roughly speaking of showing that ()0.5|) extends to Cole- 
man functions of several variables and then taking the limit at (1,1). 

Since taking the limit turned out to be rather involved in |Flj . we 
opted for an alternative approach, which was motivated by a letter 
of Deligne to the second named author D2J. Deligne observes that 
taking the limit at 1 for the multiple polylogarithm can be interpreted 
as doing analytic continuation from tangent vectors at and 1, using 
the theory of the tangential basepoint at infinity introduced in |Dlj . To 
analytically continue (jU.5|) and obtain the series shuffle product formula 
we analyze a more general notion of tangential basepoint sketched in 
loc. cit. and examine among other things its relation with Coleman 
integration. 

To give a precise meaning of the limit value to (1,1), we work over 
the moduli space A4q,5 of curves of (0, 5)-type and the normal bundles 
for the divisors at infinity -Mo,5 — -M.o,s (-^o,5 : a compactification of 
Mofi)- Two variable p-adic multiple polylogarithms are introduced. 
They are Coleman functions over *Mo,5- In their analytic contin- 
uation to the normal bundle will be discussed. In particular, we will 
relate the behavior of the analytic continuation of two variable multiple 
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polylogarithm to a normal bundle with one variable multiple polylog- 
arithm and then we get p-adic multiple zeta values as "special values" 
of two variable multiple polylogarithms. 

It is a pleasure to thank P. Deligne and A. Goncharov. The first 
named author was supported by a grant from the Israel Science founda- 
tion. The second named author is supported in part by JSPS Research 
Fellowships for Young Scientists. He would like to thank the Center 
for Advanced Studies in Mathematics at Ben-Gurion University for 
supporting a visit to Beer-Sheva in February 2003 when a large part 
of the research described in this work was carried out and he would 
also like to thank the Johns Hopkins University for their support and 
hospitality during the JAMI conference. 



I. TANNAKIAN DUALITY AND FUNDAMENTAL GROUPS 

This section is introductory, and contains some basics of the theory 
of Tannakian categories (for which the standard reference is |DMj ) and 
some ways it can be used to construct fundamental groups and path 
spaces in algebraic geometry. 

The theory of Tannakian categories has its roots in Tannaka-Krein 
duality, itself a generalization of Pontrjagin duality, which shows how 
to reconstruct a compact topological group out of the category of its 
representations. Subsequently, it was realized by Grothendieck that 
the reconstruction can be done, in the somewhat different context of 
affine group schemes, based on some formal properties of this category 
of representations. Consequently, any category with some additional 
structures and properties, formalized in the notion of a Tannakian cat- 
egory, gives rise to an affine group scheme. 

A Tannakian category (or more precisely a neutral Tannakian cat- 
egory) over a field F is a rigid abelian F-linear tensor category pos- 
sessing a fiber functor [DM] Definition 2.19]. An abelian F-linear cat- 
egory T becomes a tensor category with the addition of a bifunctor 
® : T x T — ► T satisfying certain axioms mimicking those of the ten- 
sor product of vector spaces. It is rigid if it satisfies [DM1 Definit ion 
1.7]. A tensor functor G : T — ► T is a functor together with a functo- 
rial isomorphism G(X) ® G(Y) — > G(X £g> Y) satisfying a few obvious 
properties DM Definition 1.8] 

The category VecF of finite dimensional vector spaces over F is a 
rigid abelian F-linear tensor category. The final data in the definition 
of a Tannakian category is a fiber functor, which by definition is an 
exact F-linear faithful tensor functor uj :T — > Vec^. 
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Suppose now that G is an affine group scheme over F. One than has 
the category Rep F (G) of finite dimensional F-rational representations 
of G. This is clearly a Tannakian category with fiber functor sending 
a representation to its underlying if-vector space. The main result of 
Tannakian duality [DM, Theorem 2.11] says that Rep F (G) determines 
G uniquely, and conversely, any Tannakian category arises as Rep F (G) 
for some G. 

Recall that an affine group scheme over F can be defined as a group 
valued functor on the category of F-algebras which is represented by 
an F- algebra. The way that Rep F (G) determines G is as follows: G 
is isomorphic to the group valued functor Aut®(o;) whose value on an 
F-algebra R is 

Aut®(o;)(i2) := Hom(cu ® R, u <g> R) , 

where u®R is the composition of to with the change of ring functor and 
Horn is in the category of tensor functors. We will often write 7r 1 (T, u) 
instead of Aut®(u;) to indicate the dependency on the category T. 

Tannakian duality can be used to define all kinds of "fundamental 
groups" in algebraic geometry and algebraic topology. The basic ex- 
ample to keep in mind is that of the category Locx of local system 
of finite dimensional Q-vector spaces on a connected topological space 
X. Given a choice of a point x G X, this is equivalent to the cate- 
gory of finite dimensional Q-representations of the fundamental group 
7Ti(X, x). It is clearly a Tannakian category over Q. Each point x £ X 
determines a fiber functor u x where ui x {V) = fiber of V at x. The 
obvious map iri(X,x) — > 7ri(Locx, ^i)(Q) is explicitly given by 

(1.1) 7h(1/h monodromy along the loop 7 G Aut(u x (V))) . 

We may also use Tannakian duality to define generalized path spaces. 
Namely, if T is a Tannakian category and i0\ and 102 are two fiber func- 
tors on T. Then there is a affine scheme vri(T, oj\, 0J2) = Hom^^i, u^) 
whose points in an F-algebra R are 

Horn® ,u 2 )(R) := Hom(wi <g> R,co 2 ® R) , 

Again Horn taken in the category of tensor functors. There is an obvi- 
ous composition map, 

(1.2) o : tti(T,u;2,^3) x 7Ti(T,Wi,w 2 ) -> Tri(T,u u u 3 ) , 

and in particular vri(T, ui, 002) is a homogeneous space for an action 
from the left of 7Tx(T, uji) as well as for an action from the right of 
7Ti(T, u 2 ) and is principal for both actions. For T = Locx we have a 
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map 7Ti(X, x,y) — > tti(Locx, <^x, ^y)(Q), where 717 (A, x, y) is the space 

of paths between x and y on A up to homotopy, given by 

(1.3) 

7 h. (V m. parallel transport along the path 7 G Hom^fV), k> y (V))) . 

A tensor functor F : 5 — > T clearly induces a map of path spaces 

(1.4) F* : 7Ti(T, ^1,^2) — >■ 7ri.(<S, o>i oF,w 2 o F) 

which is explicitly given by sending 7 e Horn® (cji ® -R, uo2®R) to F*{^f) 
given by F*(7)s = 7f(s) • This is clearly compatible with composition. 
In the case of categories of local systems, when F is induced by the 
continuous map / : A — > Y it is easy to see that this map is compatible, 
via ()1.3|) with the usual map induced by / on path spaces. 



2. Coleman's ^-adic integration 

In this section we recall the theory of Coleman's p-adic integration 
in several variables introduced in |Beslj and make some additional 
constructions which are required for the current work. In some respects, 
the similar theory of Vologodsky [V] is more suitable for the current 
setting. However, because it is much less explicit it seems harder to 
adapt it to the tangential constructions of £0 

Recall from |Beslj that the basic setup for Coleman integration the- 
ory is as follows: We have a field of characteristic p, k, which we 
assume for simplicity to be algebraically closed and a discrete valua- 
tion ring V with fraction field K and residue field k. A rigid triple is 
a triple T = (X, Y, P) consisting of P, a formal p-adic V-scheme, Y a 
closed K-subscheme of P which is proper over Spec(/t) and X an open 
/t-subscheme of Y such that P is smooth in a neighborhood of X. 

A simple case where rigid triples arise (see |Besl| Definition 5.1]) 



is that of tight rigid triples. Let X C y be an open immersion of 
V-schemes such that X is smooth and y is complete. The associated 
triple Trx,y) '■= (X £g> v K, y <S>y K, y), where y is the p-adic completion 
of y is called a tight rigid triple. An affine rigid triple is a tight rigid 
triple Trx,y) with X affine. 

Given a rigid triple as above we have the category of unipotent over- 
convergent isocrystals on T, denoted Un(T). This is the category of 
sheaves on the tube of X in P, overconverging into the tube of Y in 
P in the sense of Berthelot |Berj , with an integrable connection. This 
category depends, up to a unique isomorphism, only on A, so we will 
often denote it simply by Un(X). 
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A Frobenius endomorphism of X is a K-linear endomorphism of X, 
which is some power of the geometric Frobenius. A Frobenius en- 
domorphism induces an auto-functor 0* on Un(X). Suppose that 
T = T(x,y) is a tight rigid triple and that is the reduction of an en- 
domorphism ip of the pair (X, y). Then the functor 0* is given by the 
obvious functor induced by ip on Un{T). 

Definition 2.1. A Frobenius fiber functor on Un(X) is a fiber functor 
uj together with an isomorphism of fiber functors (which we write as 
an equality to simplify matters) to o 0* = uj 

For two Frobenius fiber functors U\ and u% there is an obvious action 
of 0* on the path space 7ii (JJn(X), ui, 0J2) coming from the action (jl.4j) . 
The main result of [BesTj (Corollary 3.2) can be stated, in a slightly 
generalized form, as the following: 

Theorem 2.2. Given any two Frobenius fiber functors onUn(X), uoi 
and UJ2, there exists a canonical invariant path, i.e., an isomorphism of 
fiber functors, a ul ^ W2 : UJ\ — > u 2 , which is fixed by 0*. We will call the 
isomorphism oj\ — > uj 2 analytic continuation along Frobenius. It has 
the basic compatibility property that o a WljW2 = a^^. 

Proof. We only sketch the main idea, which is rather simple, and refer 
to reader to loc. cit. for the full story. The key point is that for any 
Frobenius fiber functor u there is an action 0* on the fundamental 
group n 1 (Un(X)), which is a pro-unipotent group, and its Lie algebra 
is controlled by negative tensor powers of H} lg (X/K). The action of 0* 
on H^ ig (X/K) has strictly positive weights, and consequently one can 
show that the map 

9 >-> <t>*(9)9' 1 , ni(Un(X),uj) -> 7n(Un(X),uj) , 

is a bijection. Now, the space 7Ti(Wn(X),a>i,a>2) is a iti(Un(X), u)\)- 
principal space and its 0* action is compatible with the group action. 
From this it is very easy to deduce that there exists a unique 0*- 
invariant element. The uniqueness immediately implies the compat- 
ibility. □ 

Usually, Frobenius fiber functors as above are obtained from geo- 
metric points x : Spec k — >• X by pullback, u x = x*, assuming that 
x is fixed by 0. Concretely, for an overconvergent isocrystal (M, V) 
the fiber functor u x (M, V) is realized as the vector space of horizontal 
sections of V on the residue disc 



U x =]x[p (=points reducing to x) 



s 
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which is the tube in the sense of Berthelot. Since each point of X 
will be fixed by some Frobenius automorphism it is not hard to see 
that the invariant path is independent of the choice of the Frobenius 
endomorphism. 

An abstract Coleman function on T is defined to be a triple (M, s, h) 
where 

• M = (M, V) is a unipotent isocrystal on T. 

• s G Hom(M, O). 

• h is a collection of sections, {h x G M(U X ), x G X}, with 
V(/ij;) = 0, which correspond to each other via analytic con- 
tinuation along Frobenius. 

One can evaluate an abstract Coleman function on each residue disc 
U x by taking s(h x ). Coleman functions are then equivalence classes of 
abstract Coleman functions under a relation which essentially guaran- 
tees that their evaluations will be identical. In fact, a Coleman function 
is a connected component of the category of abstract Coleman func- 
tions, where a morphism (M\, si,h\) — ► (M2, 52,^2) is a horizontal 
map u : M\ — > M2 such that S2 ° u = s± and u(hi) = 1%2- It is shown 
Besl| Corollary 4.13] that two Coleman functions which coincide on 



m 



an open subset of the tubular neighborhood ]X[p of X in P coincide 
everywhere. 

Remark 2.3. It will be important to note how the theory can be used 
to define iterated integrals. For a similar and more detailed discus- 
sion see the proof of Theorem 4.15 in [BeslJ. Suppose that we have a 
collection of Coleman functions F^, arising from the abstract Coleman 
functions (Mi, Si, hi), and one-forms uJi for % = 1, . . . , k. Suppose now 
that the locally analytic differential form uj = Yl is closed, and 
that we wish to find a Coleman function F such that dF = u. We 
may begin by constructing the connection M = (©Mj) © 0t with the 
connection given by the formula 

V(mi, • • • , m k , f) = (Vimi, • • • , V fc m fc , df - ^ UiS^rrii)) . 

This is clearly a unipotent connection. If it is integrable we immedi- 
ately obtain our required Coleman function represented by the abstract 
Coleman function (M,s,y), where s is the projection on the last fac- 
tor and h is obtained by choosing one residue disc U, extending the 
sum of the hi to a horizontal section of M on U, and then analytically 
continuing along Frobenius. In general, V will not be integrable. How- 
ever, there is a maximal integrable subconnection (M mt , V) of (M, V) 
|Besl| (2.3)]. The fact that ^i 7 ^ is closed guarantees that we may 
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choose h to be in M mt (U) hence we obtain our Coleman function from 
(M int ,s,h). 

One can apply the theory to other fiber functors. A simple case 
was already discussed in |Beslt Section 5]. Here we will need an easy 
generalization (to the two dimensional case for simplicity) as follows: 
Suppose that T = (X, Y, P) = T(x,y) is a rigid triple and that y G 
D = Y — X is a closed point and that D is locally given near y by 
the reduction of two parameters ti and t2 on Y. Let A be the ring of 
Laurent series with coefficients in K, 



A :={/(ti, h)= 

(i,i)ez2 



/ converges on r < |ti|, 1^1 < 1 for some 



It follows from |Berj that overconvergent isocrystals on T give rise to 
a connection on a (free) A-module. 

Proposition 2.4. An overconvergent unipotent connection on T has a 
full set of solutions in B := A[log(tx), logfo)]- 

Proof. Overconvergent unipotent crystals are solved by iterated inte- 
gration and the ring B is easily seen to be closed under partial integra- 
tion with respect to ti and t 2 . □ 

Note that the ring B is independent of choice of parameters t\ and 
t 2 . The following operation will be required later. 

Definition 2.5. The constant term with respect to the parameters t\ 
and t 2 of / G B is defined as follows: Let / = fij l°g l (^i) log 7 ^)- 
Let foo = ^ciijtitff Then the constant term is aoo 

Note that the constant term does depend on the choice of parameters. 

Definition 2.6. The fiber functor uj v of Un(X) associates to a unipo- 
tent overconvergent isocrystal the vector space of solutions in B. 

Suppose that we have an endomorphism tp of the pair (X, y) reducing 
to a Frobenius endomorphism fixing y (a sufficiently high power of a 
Frobenius endomorphism will fix y). Then, pulling back by <p gives an 
isomorphism u y o ip = u y making it a Frobenius fiber functor in the 
sens of Definition 12.11 It is therefore possible to analytically continue 
along Frobenius to or from u y . 

To end this section we note that by the naturality of the analytic con- 
tinuation along Frobenius, it clearly extends to pro-unipotent isocrys- 
tals. We also note that as long as the underlying connections are defined 
over a discretely valued subfield it is possible to work over C p . 
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3. Tangential basepoints 

In this section we first recall the theory of the tangential basepoint 
in the de Rham setting due to Deligne |Dlj . In its simplest form this 
theory allows one to define fiber functors for the category of integrable 
connections on a curve associated with a tangent point "at infinity". 
We explain Deligne's interpretation of constant terms in terms of these 
tangential basepoints. In loc. cit. a higher dimensional theory is also 
sketched with no details, but these can easily be filled in, as we will do. 
We then discuss the possibility of analytically continuing solutions of a 
unipotent integrable connection along Frobenius a la Coleman to such 
fiber functors. We finally explain how one can iterate the construction 
of tangential basepoints in two different ways and get the same result 
and we close with some applications to Coleman integration theory. 

We first recall the construction of Deligne |DH 15. Theorie algebrique]. 
Let C be a curve over a field K of characteristic 0, smooth at a point 
P, and let t be a parameter at P (Deligne immediately passes to the 
completion at P but we will not do the same). Suppose (M, V) is a con- 
nection on C — {P} with logarithmic singularities along P. This means 
that locally near P the connection V can be written as V = d + T, 
where T is a section of End(M) ® fi^(logP) with at most a sim- 
ple pole at P. The parameter t induces naturally a parameter i on 
the tangent space Tp{C) (the linear parameter taking the value 1 at 
the derivation d/dt). We associate with this data the connection on 
Tp(C) — on the trivial bundle with fiber Mp (fiber of M at P), given 
by Resp(V) := d + (Resp T) dlog(t). This clearly defines a functor 
Resp. 

An easy computation shows that the functor Resp does not depend 
on the parameter t. Deligne gives an alternative, coordinate free de- 
scription of the same construction, making this fact evident. The val- 
uation vp on the fraction field K{C) gives an algebra filtration Fp 
on this field and there is an obvious canonical isomorphism between 
the associated graded algebra GipK(C) and the coordinate ring of 
T P (C) — 0, given by sending a cotangent vector dpf, thought of as a 
linear function on Tp(C), to the image of / in GrpK(C). 

Let VLQ c j K {\ogP) be the sheaf of differential forms with log singu- 
larities along P. We give &K( C y K = K{C) ®q c QQ c ^ K (\ogP) the 
filtration induced from the filtration on the first term. The nitrations 
on K(C) and on Q l K rcy K induces a filtration on M K(C) and on 
M ®0 C &K{nyK an d the assumption that V has log singularities im- 
plies that it preserves the filtration. It is easy to see that the induced 
connection on the associated graded is exactly Resp(V) 
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The construction of Resp gives us the option of producing more fiber 
functors for the category of connections on C by taking the fiber at a 
point of T P {C) — 0. A particular case of this construction gives the 
notion of a constant term for a horizontal section of this connection as 
we now explain. Suppose that the connection V is unipotent. Then it 
is very easy to see that one can find a basis of formal horizontal sections 
to V near P with coefficients in the ring K[[t]} [log(t)], where log(t) is 
treated as a formal variable whose derivative is dt/t. Let v be such a 
horizontal section. There is a sense in which we can specialize v to the 
fiber Mp, namely, taking the constant term. 

Definition 3.1. The constant term of v is obtained by formally setting 
t = log(t) = 0. 

The justification for this definition is that over the complex numbers 
one has lim^o t log(t) = 0. Thus in situations where the solutions has 
coefficients in K[[t]] + tK[[t]][\og(t)] this is indeed the constant term. 
The same can be argued p-adically, provided one takes the appropriate 
notion of limit. It is important to note that the situation described 
above is indeed what happened for p-adic multiple polylogarithms near 
1, by the main result of |F1| . 

Similarly, there is a basis of (this time global) solutions to Resp V 
with coefficients in if [log (t)]. In fact, all the solutions are of the form 
exp(Resp T ■ \ogi) ■ v with v G Mp and the exponential is a finite sum 
as RespT is nilpotent. In this case we can again take the constant 
term. This can now be interpreted as simply evaluating at 1 with the 
convention that log(l) = 0. Thus, we may formally interpret taking 
the constant term as continuing to the tangential vector 1 at P. To 
make this more than a mere heuristic, though, one needs to introduce 
a topology. It can be made precise in the complex case |Dlj and we 
will show this also in the p-adic case (see below Proposition I4.5|) . 

The higher dimensional generalization is now fairly clear. Suppose 
that a smooth variety X is given and in it a divisor D = Yli<=i Di with 
normal crossings. We assume that all the components Di are smooth. 
For a subset J C / we set Dj = HjtzjDj. Let Nj be the normal bun- 
dle to Dj and let Nj be the complement in Nj of Nji\dj for J' C J. 
Note that one only need to take J' smaller by one index and that N$ 
is considered as the zero section. Thus, for example, N® is the (one- 
dimensional) normal space to Dj minus the zero section. Finally, we 
denote by N%° the restriction of N° 7 to D°j := Dj - Uj^jDj. Follow- 
ing Deligne we construct, given a connection on X with logarithmic 
singularities along D, a connection on every Nj° with logarithmic sin- 
gularities "at infinity" . Infinity here means the union of the hyperplane 
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at infinity for the normal bundle Nj, the hyperplanes Nj — Nj, and 
Nj\ Dj - D °. 

For each j e J consider the valuation Vj on K(X) associated with 
the divisor Dj. Let 0x{D~ l ) be the localization of @x at D. There 
exists a multi-filtration Fj on 0x(-D _1 ), indexed by tuples x = (Xj)jeJ^ 
such that Fj is the 0x-module generated by {/ e 0x(-D~ 1 ), ^(/) > 
Xj for all j G J}. It is easy to see that Spec(Grj 0x(-D -1 )) is precisely 
A^j°. Again, a natural map sends a differential form df, viewed as a 
function on Nj°, linear in the bundle direction, to the image of / in 
Gy 1 j x (D- 1 ). 

Remark 3.2. Let Spec(A) be an affine patch on X with coordinates 
ti e A, i e such that fl Spec(A) is defined by ti = 0. Let B = 
A/(tj,j e J), such that Dj n Spec(A) = Spec(B) Let := Yljejt? ■ 
Then = At* and it follows easily that Grj = B[tf,j J] 

Now suppose we have a connection with logarithmic singularities 
along D, V : M -> M® 0x Q^(log J D). We give ^(D -1 ) = fi^(log£>)<8> 
0x(-D _1 ) the induced filtrations from the filtration on 0x(-D -1 )- It is 
easy to see that the differential d preserves the filtration. Now let 
M(D~ l ) = M ® QxiD- 1 ) and M ® Qx^D' 1 ) have the induced fil- 
trations. It follows that the extended connection V : M(D^ 1 ) — > 
M <S> f2^(-D _1 ) respects the filtration. The connection we are looking 
for is the gr of this connection. It is evident from this construction 
that if V is flat, so will be the resulting connection. 

Definition 3.3. We call the resulting connection the residue connec- 
tion on Nj° along Dj and denote it by Res^ j V. 

Note that the construction really depends on both D and J and not 
just on Dj. It is evident that Res^j is a tensor functor. 

As a consequence of the construction above we obtain more fiber 
functors for connections. 

Definition 3.4. If x e Nj° we let uj Xj j be the fiber functor u x>J : = 
lu x o Res DjJ . 

Suppose now that / : (X f ,D') — > (X,D) is a morphism (i.e., takes 
D' to D) and that f*Dj = UjD'-. We get an induced map 

by looking at the map on the graded sheaf of rings, obtained from /, 
sending Grj to Gr^*, where n - x is the vector with components rijXj- 
The following Lemma is trivial. 
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Lemma 3.5. Suppose that Fr is the absolute Frobenius. Then is 
also the absolute Frobenius. 

The constructions above can be iterated, and they commute in a 
sense we will now explain. For simplicity, we now assume that the 
divisor D has only two components, D\ and Consider N\. It has 
on it a divisor D' with normal crossings consisting of D[ = the zero 
section and D' 2 = restriction of N\ to D^i = D\C\D2- The complement 
of this divisor is exactly N®°. We have the following result. 

Proposition 3.6. Write N = iVS n , be the construction above per- 
formed on N\ with respect to D' . Then there exists a natural isomor- 
phism ij) : N —> N®2 in such a way that for any connection M on X 
with logarithmic singularities along D we have a natural isomorphism 
ip* ResD 5 i2 M — > ResD',12 R-es^,! M. 

Proof. The isomorphism if) is essentially the isomorphism, existing on 
any object with two nitrations, Gr 12 = Gr 2 Gr!, applied to the filtra- 
tions we have. The isomorphism on the Res's is then a completely 
formal consequence of this. □ 

Remark 3.7. To make very concrete the situation of the last Propo- 
sition, and to make it clear how it is going to be used, consider the 
simplest situation, X is the affine plane A 2 with coordinates x and y, 
D = Di U D 2 = {x = 0} U {y = 0}. Clearly, N±, the tangent space 
to {x = 0}, is again isomorphic to A 2 , and we may consider on it the 
curve C which is the section x = 1. Proposition 13.61 together with 
the functoriality of the Res construction, allows us to get the following: 
Consider the fiber functor on connections on X with logarithmic singu- 
larities along D: First apply Resi, then restrict to C, finally consider 
the fiber functor at the tangential point 1 to C at y — 0. Then this 
fiber functor is naturally isomorphic to the fiber functor at the normal 
vector (1, 1) at the point (0,0). 

4. Analytic continuation to tangential points 

We now turn to the relation between the tangential theory described 
in the previous section and the analytic continuation along Frobenius 
explained in §21 Suppose then that K is a p-adic field with ring of 
integers 0k and residue field k. Let X be a smooth scheme over &k 
and let D = Di be a divisor with relative normal crossings on X. 
We assume as before that the components Di are smooth. The schemes 
Nj° are then also smooth over 0k- Define the category Un(X, D) of 
unipotent connections on Xk with logarithmic singularities along Dk- 
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This is a Tannakian category over K. Suppose that x G (X — D)(k) and 
y G Nj°(k) are two points on the special fibers. Restricting (M, V) G 
Un(X, D) to U x and taking horizontal sections we obtain a fiber functor 
u x on Un(X, D) as in ^ Similarly, restricting Res£) Ki j V to the residue 
disk U y of y in Nj° and taking horizontal sections we obtain a fiber 
functor u y . Then we have the following result. 

Theorem 4.1. For any two points x and y, which are either in (X — 
D)(k) or in Nj°(k) for some J, there exists a canonical isomorphism 
(analytic continuation along a Frobenius invariant path), a XjV : u x ^ 
u y , such that the following properties are satisfied. 

(1) These isomorphisms are compatible in the sense that a yjZ oa Xj y = 

(2) Ifx, y G (X—D)(k) then a XiV is the isomorphism of Theorem ] 2. 21 
pulled back via Un(X, D) — > Un((X — D) K ) 

(3) Ifx, y G Nj (k>) (for the same J) then a x<y is the isomorphism of 
Theorem\EEpulled back via Res Dif ,j : Un(X, D) -> Un((N?) H ) 

The Frobenius invariant path is functorial with respect to morphisms 
(X', D') — > (X, D) in the obvious sense. 

Proof. Since the nerve of an affine covering is contractible, we are easily 
reduced to the following case: X = Spec (A) is affine and D is defined 
on it by the vanishing of t\ • t% ■ ■ ■ t^, where ti are parameters. Let A 1 * 
be the weak completion of A in the sense of [Put and let be the 
weak completion of A[(t\ ■ ■ -tfc) -1 ]. We clearly have A* C B\ We may 
further assume that there exists a map, whose reduction is a Frobenius 
endomorphism, (p : B' — > B* preserving A* and sending ti to t\. We 
consider the de Rham complex of A^ with logarithmic singularities 
along D K , 

VTZ(A\ D) : Aj, ± ^^(Ak/AOdlogfo) i> ■ ■ ■ 

and the de Rham complex VTZ^B^) of B^ K . Note that the cohomology 
of this last complex is the Monsky-Washnitzer = rigid cohomology 
H rig ((X — D) K ). There is an obvious map 

(4.1) VK{A^D)^VK{B\ { ) 

and an action of tp* on both complexes compatible with this map. 

Proposition 4.2. The map ()4.1|) is a quasi-isomorphism. 

Proof. This follows essentially from the arguments in |BC| Section 6] 

□ 
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We only need the fact that the induced map on H 1 is an injection. 
This immediately implies the following. 

Corollary 4.3. The map if* : H l (VK{A^D)) -> H\VH{A^,D)) is 
independent of the choice of if and has strictly positive weights, 

Proof. The same is true for H^ ig ((X — D) K ) as in the proof of Theo- 
rem o 

□ 

Define now the category Un(A>, D) of unipotent integrable ^-connections 
with logarithmic singularities along D. There is an action of tp* on this 
category, hence on its fundamental group and torsors. By the same 
argument as in the proof of Theorem 12.21 we have using Corollary 14.31 
that there exists a unique invariant path between any two fiber func- 
tors on Un(A',D) and that this is independent of the choice of (p. 
We have such a fiber functor uo x for x G (X — D) K , and since there 
is an obvious residual functor Un(A\D) — > Un((Nj ) K ), defined in 
the same way as in Definition 13.31 we also have fiber functors u y for 
y G (Nj°) K ). The fiber functors discussed before the theorem are the 
compositions of the ones we just defined with the restriction functor 
Un(X,D) — > Un(A\D) and we obtain our required path by pulling 
back along this restriction. With this definition, and (J2J) are essen- 
tially clear, while (jSJ) follows easily because by Lemma 1331 we have that 
if induces residually a morphism on Nj° whose reduction is a Frobenius 
automorphism, and consequently the pullback via the residue functor 
of a Frobenius invariant path is (^-invariant. □ 

We can use the extension of Coleman functions described here to 
make a theory of Coleman functions "of algebraic origin" on the pair 
(X, D). The abstract Coleman functions are then triples (M, s, h) with 
M G Un(X, D) and s : M — > 0x K , while h is a compatible system of 
local horizontal sections as before. Theorem 14.11 allows us to define a 
map Res£>,j from these algebraic Coleman functions to Coleman func- 
tions on the normal bundle Nj°. It is defined by 

Res AJ (M, s, h) = (Res DJ M, Gr s, h) 

where h now refers to the induced compatible system of horizontal 
sections whose existence is guaranteed by Theorem 14.11 

In particular, we obtain from a Coleman function / as above a Cole- 
man function f^ on the normal bundle to D. The following proposi- 
tions suggests how to compute it. 

Proposition 4.4. Suppose that we have the relation df = Y^^iQi with 
f and gi Coleman functions of algebraic origin. Then we have the 
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relation df^ = ^(Resp, uji)g\ D \ where, if uj is locally written as uj' + 
/idlog(t), with t the defining parameter for D, then Res£>(cG>) = oj'\d + 
h\ D d\og(t). 

Proof. We consider the way / can be written as a Coleman function 
following Remark |2. HI In that description suppose that the connection 
M constructed there was indeed integrable. Then the result is easily 
obtained by computing the residual connection for M. The result re- 
mains true if M is not integrable. To see this we first note that the M mt 
construction is algebraic by Lemma 2.4 in |Bes2j and the discussion fol- 



lowing it. One further sees that it is a connection with logarithmic sin- 
gularities along Dk- It follows that we may obtain / from the abstract 
Coleman function (M mt , s, h) as in Remark 12.31 The residual construc- 
tion is not limited to integrable connections and is functorial. Thus 
there is a horizontal map of connections ResD,i(M mt ) — > Res^^M). 
Consequently, f^ D \ which is defined on the normal bundle to Di by 
Res d, i(M mt ,Gr s, h), still satisfies the differential equation given by 
Res D /(M) □ 

We can now give, in the p-adic tangential basepoint interpre- 

tation of the constant term. 

Proposition 4.5. In the situation of Definition \2. 5\ the constant term 
with respect to the parameters ti and t 2 of a horizontal section of a 
unipotent overconvergent isocrystal V on T is the evaluation at the 
normal vector t\ — 1,?2 — 1 °f the analytic continuation along Frobe- 
nius of this horizontal section. 

Proof. We only prove the 1-dimensional analogue as the proof is similar. 
In fact, we prove something stronger. Suppose (using the notation of 
this section) that the connection is locally given by d + V. Let / be a 
horizontal section for the connection V and let v be its constant term. 
We claim that the analytic continuation of / is precisely the solution 
c(f) := exp(RespT ■ logf) ■ v (the result is then proved by specializing 
to t — 1). We already know that this is indeed a solution of Resp(V). 
The association V i— > (/ i— ► c(/)) is clearly compatible with direct 
sums and tensor products hence defines a path. To show that this is a 
Frobenius invariant path we may assume that ip is such that <p*(t) = t p 
The constant term of <£>*(/) remains v while y?*V is associated with 
ip*T. It follows that 

c(<p*f) = exp(Res P (^T)-logf)-f = exp{pRes P T -logt) ■ v = c(f)(i p ) , 
which proves our claim. □ 
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5. The set up of the moduli space Mo,5 

This section is devoted to giving a quick review on basic materials 
of the moduli space M.0,5 of genus curves with 5 distinct marked 
points and its Deligne-Mumford compactification Mo,5- The two (and 
one) variable p-adic multiple polylogarithms are introduced as Coleman 
functions on .Mo 5- 

The moduli space M , 5 = {(P,)f =1 G (P 1 ) 5 !^ + Pj (i + j)}/PGL(2) 
is identified with 



(5.1) {(x,y) G A 2 }\{x = 0}U{y = 0}U{x = l}U{y = l}U{xy = 1}. 

This identification is given by sending (x, y) to 5 marked points in 
P 1 given by (0, x, 1, ^, 00). The symmetric group S 5 acts on A4 0t5 by 
a(P i ) = P a -i {i) (1 < % ^ 5) for a G S 5 . Especially for c = (1, 3, 5, 2, 4) G 
S 5 its action is described by x 1— > jE^: V *~ ¥ x - 

The Deligne-Mumford compactification of Aio,5 is denoted by Aio,5- 
This space classifies stable curves of (0, 5)-type and the above ^-action 
extends to the action on A4 0t5 . This space is the blow-up of (P*) 2 (d 
^0,5) at (x, y) = (1,1), (0,oo) and (00, 0). The complement .Mo,5 — 
-M ,5 is a divisor with 10 components: {x = 0}, {y = 0}, {x = 1}, {y = 
1}, {xy = 1}, {x = 00}, {y = 00} and 3 exceptional divisors obtained 
by blowing up at (x, y) = (1,1), (00, 0) and (0, 00). In particular 
for our convenience we denote {y = 0}, {x = 1}, the exceptional 
divisor at (1,1), {y = 1} and {x = 0} by D\, D 2 , D 3 , D A and D 5 
(or sometimes D ) respectively. It is because c l (D ) = Di. These five 
divisors form a pentagon and we denote each vertex Di n D^\ by P«. 
Hence we have c 1 (Pq) = Pi. The two dimensional affine space U\ = 
SpecQ[x, y] gives an open affine subset of -Mo,5- The S'5-action gives 
other open subsets Ui = c l_1 (L r i) = SpecQ[zi,Wi] (1 ^ % ^ 5) in -M ,5 
where (z 1 ,w 1 ) = (x,y), (z 2 ,w 2 ) = (y,jE^), (z 3 ,w 3 ) = - xy), 

(z 4 , w 4 ) = (1 - xy, ij) and (z 5 , w 5 ) = x). 

For a = (a 1 , ■ ■ ■ , a k ) G Z k >Q , b = (b 1 , ■ ■ ■ , b t ) G Z l >0 , and x, y G Q p 
with \x\ p < 1 and \y\ p < 1 we define two variable p-adic multiple 
polylogarithm by 



x mfe y ni 

Lia, b (a:,J/):= \^ — ^ ^eQ p [[x,y}}, 

<ni<---<n ; 



18 AMNON BESSER AND HIDEKAZU FURUSHO 

and for c = (ci, • • • , qJ G Z> one variable p-adic multiple poly- 
logarithm by 

yrn h 

Li c (y):= £ y gQ P [[y]]cQ^,y]]. 

// &i '''llli 

0<mi<-<m k 1 ft 

It is easy to see that these functions satisfy the following differential 
equations. 

(5.2) 



d 

dx 



Li aib (x,y) 



Li ( oi,-,o fc _ 1 ,o fc -i),b(a:,J/) if Ofc 7^ 1, 

lZ^Ll( aii ... )afc _ l))b (x,y) - (i + Ll(oi,-,a fc _i,6i),(6 2 ,-,&()( X '?/) 

if a fc = 1, k ^ 1,1 ^ 1, 
i^Lib(y) - (s + T3^) Li(6i) ) (& 2) - ) 6,)( a: >S') if a fc = 1, A; = 1,/ ^ 1, 

if a fc = l,ife ^ 1,Z = 1, 
[ ^ Li b (y) - (I + ^) Li b (xy) if a k = 1, k = 1, 1 = 1, 

lU gLj{blt ... M _ ubl _ 1) (x,y) if k ^ 1, 



1 T : 

l-X ( a l:-" > a fc- 



— Li a>b (x,j/) = < J L ^U a){bl ... M _ l) (x,y) if 6, = 1,1^ 1, 
jzzU^xy) if 6/ = 1,1 = 1, 



(5.3) 



Li c (:n/) 



-U c (xy) 



(5.4) A Lic ( y ) = 




.i,c h -i)(a^/) ifc ft ^l, 



(ci,— ,Ch 



Li (ci,-,c h _i)(a^) if c fe = 1, /i ^ 1, 
if c h = 1, /i = 1, 

=i,-,c h _i,c h -i)(»y) ifa^i, 

Li (ci,-,c h _i)(a^/) ^c h = l,h^l, 
iic h = l,h= 1, 



lU {cu ... jCh _ uCh _ 1) (y) iic h ^ 1, 

T=^ Li (ci 1 - 1 c h _ 1 )(j/) ifc h = l,/i^l, 

if c/j = 1, /i = 1, 



i 

1-2/ 
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By the above differential equations, Li^^x, y), Li c (xy) and Li c (y) 
are all iterated integrals of ^, & and differential 
forms over A^o,5- Whence they are obtained from some triple (M, s, y) 
over -Mo,5- We interpret them as Coleman functions over the rigid 
triple (^0,5,^0,5)- This means that they are analytically continued 
to -Mo,5(Qp) as Coleman functions by the methods of analytically con- 
tinuation along Frobenius in £J2J 



6. Analytic continuation of two variable ^-adic multiple 

polylogarithms 

In this section two (and one) variable p-adic multiple polylogarithms 
are analytically continued into iV™ (i G Z/5), which is a Zariski open 
subset of the normal bundle N^- of the divisor Di introduced in the 
previous section. 

Notation 6.1. For a Coleman function / over M.0.5, f^ Di ^ means the 
analytic continuation of / to (i G Z/5). For a = (a 1; • ■ ■ , a^) G Z> 
and b = (b%, ■ ■ ■ ,bi) G Z' >0 , F ah stands for the Coleman function 
Li a ,b(%, y) — Li a b(xy) and for c = (ci, • • • , qJ G Z> , l7 c stands for the 
Coleman function Li c (xy) — Li c (y) over .Mo,5- 

Lemma 6.2. F^^ = and Gc = for any index a ; b and c. 

Proof . The constant terms of Li a ^(x, y), Li c (xy) and Li c (y) at the 
origin P5 are zero because there are no constant terms in their power 
series expansions. We take their differentials and take their residues at 
y = 0. It gives by induction because each term will be a multiple 
poly logarithm with one lower weight than the original one. Whence 
Li at h{x,y), Li c (xy) and Li c (y) are identically zero. It gives our claim. 

□ 

Lemma 6.3. F^^ = if a is admissible 1 and Gi? 2 ^ = for any 
index c. 

Proof . On the affine coordinate (z 2 ,w 2 ) for U 2 , the divisor D 2 is 
defined by w 2 = 0- By using dx = ^^—^Xdz 2 + , Z2 ~\^ 2 dw 2 and 
dy = dz 2 , we obtain the following from the differential equations in 



1 An index a = (ai, • • • , a/t) (aj € N) is called admissible if a% > 1- 
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(GE2D ~ (Q. 



A 

dz' 



Li a>h {x,y) = 



dz 2 



w 2 



1 — 2 2 lt>2 



^i(ai,-,a Jk _i,a*-l),b(a?,2/) + T 2 U *^,~ .ft^-l) 0, 2/) 

if a fc > 1,6; ^ 1, 



1 -^VCbi,-^-!)^'?/) 

if CL k 



if a fc > l,bi = 1,1 ^ I, 

IZ ^LZ(a 1 ,...,a fc _ 1 ,a fc -l),b(z,y) + JZ^Li^Xy) 

if a fc > 1,6/ = 1,Z = 1, 



V 

T=f^ + i) ■ Li (c!,- ,c fc -i,c h -i) M if c h ^ 1, 

(22^2-1) Z2 — 1 / 



^2"^2yj-~^2; 1 W2 — 1 

^ (Z2— 1)(Z2W2 — 1) 22 — 1 

22W 2(1— 102) 1 m 2 — 1 

)(Z2W2 — 1) Z2 — 1 



(ci,- ,Cft- 



(«2 — 1)(«2«2— 1) 22 



x )(rrj/) if c fc = 1, /i / 1, 
if c h = 1, h = 1, 



22 

T^Z 2 

1-22 



if c fe = 1, /i = 1. 



Following Proposition 14.41 we compute the residue around D 2 



d 



^Li?>\xy) 
dz 2 



-Li 



(D 2 ) 



I-22 a,(6i,' 



1 



1-22 

22 (ci,— ,c h _i,c h -l) 

1-22 [Cl,— ,C f , 



if a fc > 1, k 
if a fc > 1, bi 



1,1^1, 
1, 



1 



1-22 



(xj/) if c h ^ 1, 
^(zy) if c h = l,h=£l, 
if ^ 



1,/i 



(Cl,— ,C h _l,Cfc-l)w / ™ 7 ' 



22 



1-22 



if c h = l,h^ 1, 
if c ft = l,/i = 1. 
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Similarly we compute 
d 



dw 2 



Li^\x,y)=Q if a fc >l, 



d -Lif*\xy) = 0, 



dlU2 



d -Li D *\y) = o. 



dlU2 



These computations imply that differentials of F^ 2 ^ and Gc, with 
respect to z~2 and W2 are zero by induction. Therefore they must be 
constant. By Lemma 16.21 their constant terms at Pi is zero. So they 
are identically zero. □ 

Lemma 6.4. F^^ = if a and b are admissible and G^ 3 ^ = if c 
is admissible. 

Proof . On the affine coordinate (2:3,103) for U3, the divisor D 3 is de- 
fined by w 3 = 0. By using dx = -w 3 dz 3 -z 3 dw 3 and dy = j^^dz 3 + 
(23ffi3~-i) 2 ^ 3 ' we obtain the following from the differential equations in 

— Li^ h {x,y) = — ^-^rI«(tti,...,o jb _ 1 ,aj k -i),b(«,!/) 
az 3 z 3 w 3 — 1 

+ 1 Li a! ( 6l) ... i6i _ li6! _i)(a;,2/) if a k > 1 and 6/ > 1, 

1 — z 3 w 3 

-^-Li c {xy) =0 if c h > 1, 

17 Lic ^ = 1 ~ ^(ci,-,^-!,^-!)^) if c ft > 1. 

Following Proposition \4A\ we compute the residue around D 3 

-^-L^ D b 3) (x, y) = if a fe > 1 and 6; > 1, 

az 3 



-^-Li^Xxy) = if c/i > 1, 

"^3 

-^—Li^P 3 \y) = if c ft >l. 

etz 3 
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Similarly we also compute 

-Li^\x,y) = if a& > 1 and bi > 1, 



d,W3 



d -Lif*\xy) = ifc h >l, 



dlU3 



d -LiV } *\y) = iic h >1. 



G?U> 3 



These computations imply that differentials of and Gc^' 3 '' with 

respect to z% and W3 are zero by induction. Therefore they must be 
constant. By Lemma f6. 31 their constant term at P2 is zero. So they are 
identically zero. □ 

In (Flj it was shown that the limit (in a certain way) to z = 1 
of Lijt lj ...,fc m (^), which is a Coleman function over P 1 \{0, 1, 00}, exists 
when k m > 1 (loc. cit. Theorem 2.18) and p-adic multiple zeta value 
C P (ki, • • - , k m ) is defined to be this limit value (loc. cit. Definition 2.17), 
but by using the terminologies in we reformulate its definition as 
follows. 

Definition 6.5. For k m > 1, thep-adic multiple zeta value ( p (ki, ■ ■ ■ ,k. 
is the constant term of ^ki,- ,k m ( z ) aX z = 1. 



In the case for k m = 1, the constant term of Lh^... ,k m (z) at z = 1 is 
actually equal to the (canonical) regularization (—l) m I p (BA km - 1 ~ 1 B ■ ■ ■ A kl 
of p-adic multiple zeta values by loc. cit. Theorem 2.22 (for this nota- 
tion, see loc. cit. Theorem 3.30). 

The following is important to prove double shuffle relations for p-adic 
multiple zeta values. 

Proposition 6.6. (1) The analytic continuation Li^\x,y) is con- 
stant and equal to £ p (a, b) if a and h are admissible. 
(2) The analytic continuation Li^P 3 \xy) andLi[. D3 \y) are constant 
and take value Cp( c ) if c is admissible. 

Proof . By Lemma 16.41 it is enough to prove this for Lif A \y). By 
the argument in Lemma 16.21 L%j (y) = 0. By the computation in 
Lemma 16.31 Lii D2 ^(v) = Lii D2 \z~2). So the constant term of Lic\y) 
at P2 is equal to the constant term of Li c (z2) at z~2 = 1, which is 
C P (c). By the computation in Lemma IB~4l Lii ^ (v) must be constant 
if Ch > 1. Since this constant term must be the constant term of 
Li^\y), U D3 \y) = ( p (c) for c h >l. □ 
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By discussing on the opposite divisors D 5 , D 4 and D 3 , we also obtain 
the following. 

Proposition 6.7. (1) The analytic continuation Li^\y,x) is con- 
stant and equal to £ p (a, b) if a. and b are admissible. 
(2) The analytic continuation Li[P 3 \xy) andLi^P 3 \x) are constant 
and equal to ( p {c) if c is admissible. 

7. The double shuffle relations 

In this section, we prove double shuffle relations for p-adic multi- 
ple zeta values ( Definition lfi.5j) . Firstly we recall double shuffle re- 
lations for complex multiple zeta values. Let a = (ai, • • • , a&) and 
b = (pi, ■ ■ ■ , bi) be admissible indices (i.e. > 1 and bi > 1). The se- 
ries shuffle product formulas (called by harmonic product formulas 
in |Flj and first shuffle relations in |Glj ) are relations 

(7-1) C(a)-C(b)= CWa,b)) 

which is obtained by expanding the summation on the left hand side 
into the summation which give multiple zeta values. Here 



Sh<(k,l) := \j{o- :{!,■■■ ,k + l}^ {I,-- - ,N} 

N 



a is onto, 



o-(l) < •■■ < a(k),cr{k + l) < ■■■ < <r(k + l)\ 



and cr(a, b) = (ci, • • • , cn) where N is the cardinality of the image of 
a and 

a s + b t -k if = { s it] with s < t, 

a s if (T _1 (i) = {s} with s ^ k. 

b s _ k if cr _1 (z) = {s} with s > k. 

One of the easiest example of (j7.1J) is (|0.2jl . 

On the other hand, multiple zeta values admit an iterated integral 
expression (cf. |Glj . |IKZj see also |F0j ) 



Z 1 <iu dw c?u 
(,(a) = / — o • ■ • o — o o — o o 



n u u 1 — u u 1 — u 



o — o • ■ ■ o — o 



u u 1 — u 

-v ' 

"1 



24 



AMNON BESSER AND HIDEKAZU FURUSHO 



Here for differential 1-forms ui,U2, ■ ■ ■ ,uj n (n ^ 1) on C an iterated 
integral J Q U\ o lu 2 ° ■ • ■ ° uj n is defined inductively as o>i(ti) J^ 1 co>2 ° 
• • -ou n . There are the well-known shuffle product formulas (for example 
see loc. cit.) of iterated integration 

/ Ul O ■ ■ - O UJ k • / LU k+1 O ■ ■ ■ O U k+ l — / CJ T (l) O • ■ ■ O U T ( k+ l), 

J ° J ° T£Sh(k,l) J ° 

where Sh(k, I) is the set of shuffles defined by 

Sh(k, I) :— |r :{1, • • • , k + 1} — > {1, • • • , k + 1} r is bijective, 

r(l) < • •■ < T(k),r{k + 1) < ■ ■ ■ < r(fc + Z)}. 

They induce the iterated integral shuffle produce formulas (called 
by shuffle product formulas simply in |Flj and second shuffle relations 
in for multiple zeta values 

(7.2) C(a)-C(b)= E CC/rO^wu) 

TeSh(N a ,N b ) 

where N a = ax~\ h a k , N h = bi H h fy. For c = (ci, • ■ ■ , c&) with 

h, ci, . . . , Ch ^ 1 the symbol means a word A Ch ~ 1 BA Ch ~ 1 ~ 1 B ■ ■ ■ A Cl ^B 
and conversely for given such W we denote its corresponding index 
by Iw- For words, W = X\---X k and W = X k+ \ ■ • ■ X k +i with 
Xi E {A, B}, and r e Sh(k,l) the symbol r(W / ", W 7 ') stands for the 
word Zi • • ■ with = J T -i(j). One of the easiest example of (|7.2j) 

is (B - 

The double shuffle relations for multiple zeta values are linear 
relations which are obtained by combining two shuffle relations (|7.3j) . 

i. e. series shuffle product formulas (j7.1j) and iterated integral shuffle 
produce formulas ()7.2|) 

(7.3) £ CHa,b))= £ C(^ aiWb ))- 

aGSh<(k,l) T£Sh{N R ,N b ) 

The following is the easiest example of the double shuffle relations 
obtained from @T2J) and (JQJ): 

((knk 2 ) + ({k2,k 1 ) + ((k 1 + k 2 ) 

i=o ^ ? ' i=o \ 3 / 

for fci, fc 2 > 1- 
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Theorem 7.1. p-adic multiple zeta values in convergent case (i.e. for 
admissible indices) satisfy the series shuffle product formulas, i.e. 

(7.4) C P (a)-C P (b)= Yl C P Ka,b)) 

creSh^(k,l) 

for admissible indices a and b. 

Proof . Put a = (oi, • • • , a*.) and b = (&i, - • • ,bi). By the power 
series expansion of Li a ^(x, y) and Li a (x) in we obtain the following 
formula 

(7.5) Li a (x) • Li b (j/) = Li a,b(^)- 

aeSh^(k,l) 

Here 

(mi,— ,m k ,ni,--- ,ni)GZ^_ + 1 K i i 

with 

^++ = {(ci, • • • , c fc+J ) e Z*J J | Cj < 9 if o-(z) < cr(j), Cj = Cj if = a(j)} 

Then for each a e Sh^(k,l), Li^ b (x, y) can be written Li a ' jb /(x, y), 
Li a ',b'(z/? x ) or Li a ',b'(^y) for some indices a' and b'. We note that, if 
a and b are admissible, then these a' and b' are also admissible. By 
Proposition 16.61 and Proposition 16.71 we know that analytic continu- 
ations Liif(x,y), 13™ (y, x), U™(xy), U™(x) and Ii£*>(y) are 
all constant and take values ( p (&, b), ( p (h, a), Cp( a ? b), Cp( a ) an d CpO 3 ) 
respectively when a and b are admissible. Therefore by taking an an- 
alytic continuation along Frobenius of both hands sides of (J7.5)) into 
iV^,(Qp), we obtain the series shuffle product formulas (|7.4j) for p-adic 
multiple zeta value in convergent case. □ 

By this theorem we say for example 

( P {h) ■ ( P {fa) = C P {h, fa) + C P (fa, fa) + C P {fa + fa) 
for fa,fa> 1 which is a p-adic analogue of (|0.2j) . 

Corollary 7.2. p-adic multiple zeta values in convergent case satisfy 
double shuffle relations. Namely 

5^ e P 0(a,b))= Y C P (^r(w a ,iy b ))- 

a£Sh^(k,l) TeS/i(JV a ,W b ) 

holds for a*, > 1 and bi > 1 . 
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Proof . It was shown in |Flj Corollary 3.46 that p-adic multiple zeta 
values satisfy iterated integral shuffle product formulas 

(7.6) C P (a)-C P (b)= Yl Cp(-fr(w.,Wb))- 

re5ft(JV a> iV b ) 

By combining it with Theorem 17.11 we obtain double shuffle relations 
for p-adic multiple zeta values. □ 

Therefore we say for example 
( p (ki, k 2 ) + ( p (k 2 , ki) + C P (h + k 2 ) 

i=o \ z ' j=o \ 3 / 

for ki,k 2 > 1 which is a p-adic analogue of (|U.4|) . 

In complex case there are two regularizations of multiple zeta values 
in divergent case, integral regularization and power series regularization 
(sec [IKZ], |Glj $2.9 and §2.10). The first ones satisfy iterated integral 
shuffle product formulas, the second ones satisfy series shuffle product 
formulas and these two regularizations are related by regularization re- 
lations. Actually these provide new type of relations among multiple 
zeta values. In the case of p-adic multiple zeta values, p-adic analogue 
of integral regularization appear on coefficients of p-adic Drinfel'd as- 
sociator (see |Flj ) and they satisfy iterated integral shuffle product 
formulas like ()7.6|) . On the other hand, it is not clear at all to say 
that p-adic analogue of power series regularization satisfy series shuffle 
product formulas and regularization relation. It is because that in the 
complex case the definition of this regularization and the proof of their 
series shuffle product formulas and regularization relation essentially 
based on the asymptotic behaviors of power series summations of mul- 
tiple zeta values (see |Glj Proposition 2.19) however in the p-adic case 
our p-adic multiple zeta values do not have power series sum expression 
like (jU.lj) . Recently the validity of these type of relations among p-adic 
multiple zeta values were achieved in |FJj by using several variable p- 
adic multiple polylogarithm and a stratification of the moduli A4q n+3 
(N^3). 
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